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FOREWORD 


V 


This document describes the processing software programs developed under the Marine 
Seismic Display Enhancement Program. These programs provide the geophysicists and 
acousticians with software to enhance marine seismic profiling data. This report contains 
a general software program description, UNIVAC 1108 computer listing, and the operation 
procedure for each program. Software to decode the Seismic Data Acquisition System's 
data, to process the data, and to prepare the data for playback on the Seismic Data Display 
System is included. Required modifications to the UNIVAC 1108 computer listings for 
operation on the CDC 6600 are given. “Ben oa Pe che) 


This document is Volume II of the Marine Seismic Display Enhancement Program Final 
Report. - ‘Volume I presents a summary of project results with examples of processed 
seismic profiling data. Volume III describes the operation of the hardware acquisition and 
display systems developed under this program. 
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Loa eer ABSTRACT 


This dooumentexplains six software programs developed to enhance display of digitally 
recorded seismic data. The first two programs are digital filters which remove undesired 
frequency content from seismic data, he third vertically stacks several records 
with the same source to hydrophone array distance resulting in improved signal-to-noise. 
The fourth program is an optimum least squares filter commonly known as the t wiener 
Filter# 3 he fifth § gevelops a band limited zero-phase time domain pulse with 
defined frequency content;? Phe sixth progranf develops a spectral analysis routine for time 
series data by use of a fast fourier transform method. complete listing of the program 
for implementation on the UNIVAC 1108 is provided andthe required changes for implemen- 
tation on the CDC 6600 are given. sas 
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This document is the result of the Marine Seismic Enhancement Program (1977). 
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MARINE SEISMIC DISPLAY ENHANCEMENT PROGRAM 
PROCESSING SOFTWARE 


I. INTRODUCTION 


This document describes the software programs developed to process seismic pro- 
filing data. These software programs were developed because seismic systems used to 
gather data introduce data distortions. These distortions may be frequency nulls, phase 
shifts, polarity reversals, secondary pulses, pulse images or multiples. Each of these 
distortions, when introduced into the data, reduce resolution and signal-to-noise. These 
programs are designed to remove these data distortions, thereby improving data quality. 


This is Volume II of a three series set. Volume I presents a summary of project 
results with examples of processed and displayed data. Volume III explains the operation of 
the Seismic Data Acquisition System (SDAS) and the Seismic Data Display System (SDDS). 


Six main programs have been implemented for processing seismic profiling data 
digitized by the SDAS (Seismic Data Acquisition System). These programs are: 


TITLE PROGRAM NAME 
1) Time Domain Filter MAINO1 
2) Frequency Domain Filter MAIN02 
3) Vertical Record Stacking Program MAIN03 
4) Wiener Filter Program MAINO4 
5) Desired Waveshape Program MAINO7 
6) Spectral Analysis Program SDAS 


The programs read an SDAS formatted tape having a maximum of 6000 data points per 
record. These programs operate on two computer systems; the UNIVAC 1108 or the CDC 
6600 interfaced thru a CDC 1700 computer. 


I. PROGRAM DESCRIPTION 


A. TIME DOMAIN FILTER (MAINO1) 


The Time Domain Filter (MAINO1) uses a zero phase Linnette Filter (Linnette, 1961). 
There are four filtering options: 


1) Lowpass Filter: Passes all frequencies below a cut-off frequency 
2) Highpass Filter: Passes all frequencies above a cut-off frequency 


3) Bandpass Filter: Passes all frequencies between a high cut-off frequency and 
a low cut-off frequency 


4) No Filter: Passes all frequencies unattenuated 
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This program uses three input parameters for the Highpass Filter or the Lowpass 
Filter. These parameters are CUT, H, and N. CUT relates to the cut-off frequency (Fo) of 
the filter. H is the slope of the filter. The actual cut-off frequency of the filter is 
Fo = CUT + H. N indicates one-half of the number of filter weights used. A maximum of 
300 filter weights is permissible. 


Bandpass filtered data is produced by highpass filtering the data, then lowpass filter- 
ing the data. Filter parameters are entered separately for the Lowpass filter and the 
Highpass filter. 


The option of ''no filter'' produces an unprocessed display of the data. 


B. FREQUENCY DOMAIN FILTER (MAINO2) 


This filter converts the data to the frequency domain employing a fast fourier trans- 
form. The program then multiplies this data by a tapered step function (controlled by the 
bandwidth). The available tapered step functions are: 


1) rectangular (no taper) 

2) cosine 

3) gaussian 

4) triangular (linear taper with frequency) 
The resulting data is converted back to the time domain via an inverse fast fourier transform 
resulting in bandpass filtered data. 


This filter is designed to bandpass filter data but it can be forced to lowpass filter 
data by setting the input parameter FRQLOW equal to the negative of the parameter FRQHI. 


At present the program will filter a maximum of 2048 data points. To filter larger 
data records the following changes must be made. : 

1) dimension of ''Y"' in MAINO1 increased 

2) dimension of ''P' in MAINO1 increased 

3) Do loop in SUB012 Line 9 must be increased 
The number used in the dimension statements must be 2N where N is the number used in 
change 3. 
Cc. VERTICAL RECORD STACKING PROGRAM (MAINO03) 


This program stacks consecutive records. It is not meant to be used for multiple 
channel streamer data as it makes no attempt to correct for the move-out condition present 
in multi-channel streamer data. This program reads consecutive records averaging point 
by point the corresponding time period samples. 


After the number of records read equals the number of records to be stacked, output 
commences. Each new record read replaces the oldest record in the stacking process and 
the new record's "header data” is used as the "header data'' for the output record. 
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The maximum output record length is limited by the number of arrays stacked. The 
maximum output record length equals 20010 data points divided by the number of arrays. If 
stacking 10 records, the maximum number of output data points allowed is 2001. To increase 
this limitation, two values must be changed in SUB029. The first is the dimension of DDATA 
and the second is the number 20010 in line 9. These two numbers must always be equal. 


D. | WIENER FILTER PROGRAM (MAIN04) 


The Wiener Filter used consists of an input signal, a desired output signal, and an 
actual output series (Robinson, 1967). The process minimizes the power existing between 
the desired and actual output signals, resulting in the least squares optimum filter. This 
program uses two input waveshapes, the BASIC WAVESHAPE and the DESIRED WAVESHAPE, 
to develop a set of filter weights. The BASIC WAVESHAPE is a representation of the propa- 
gating sound source wavelet. The BASIC WAVESHAPE is stored in the array B and has LB 
number of data points. The DESIRED WAVESHAPE, stored in the array D, is a bandpass, 
zero phase, pulse. Its spectral content should approach that of the BASIC WAVESHAPE. 

LD is the number of points in the DESIRED WAVESHAPE. 


The input data, BASIC WAVESHAPE and DESIRED WAVESHAPE, can be manipulated 
in the following ways before developing the filter weights: 
1) BASIC WAVESHAPE may be read from any SDAS data tape. 
a) Consecutive files can be stacked to develop BASIC WAVESHAPE. 
b) The program can select the first data point of the basic from the data file. 
c) The operator can select the first data point of the basic from the data file. 


d) The period of the bubble pulse (air gun source) can be lengthened or 
shortened while maintaining the frequency content of the main pulse 
and the bubble pulse (total frequency spectrum will change because 
bubble period changes). 


2) BASIC WAVESHAPE can be read from card images (a record or records still 
read from tape). 


3) The DESIRED WAVESHAPE can be time shifted in reference to BASIC 
WAVESHAPE. 


The plot output format has two sections. The first section contains the input parameters 
used in developing the wiener filter weights whereas the second section is the resulting 
Wiener Filter processed data. 

The first section has the following information plotted: 

1) Each record used in the stack (even if basic is read in from card images) 

2) Results of stacking records (if stack greater than one) 

3) DESIRED WAVESHAPE 

4) BASIC WAVESHAPE 

5) Convolution of BASIC WAVESHAPE with" Filter Weights" 


6) Filter Weights 
7) Convolution of " Filter Weights" with stacked file 


The second section is a series of seismic records consisting of one record for each 
record processed by the filter weights, developed in the first section. 


E. DESIRED WAVESHAPE PROGRAM (MAINO7) 


This program develops a zero phase, band limited pulse waveform. This waveform 
is developed in the frequency domain and converted to the time domain. One hundred data 
points are generated to represent the DESIRED WAVESHAPE. This program writes the 
DESIRED WAVESHAPE in three locations: 


1) Memory file 15 of computer 
2) On computer listing 
3) On card output (if requested) 


This program also writes, on the computer listing, the spectral amplitudes (db) of the 
DESIRED WAVESHAPE. The bandwidth (Hz) of each Frequency bin is equal the sample 
Frequency divided by 2N, 


F. SPECTRAL ANALYSIS PROGRAM (SDADS) 


This program reads a given data record and computes its frequency spectrum. The 
frequency spectrum data is plotted by CALCOMP for presentation. The data record can be 
data directly from an SDAS tape, computer filtered SDAS tape, or card images of data. The 
program will initiate the frequency analysis on either a first return (water-sediment inter- 
face) or prechosen point within the record. The program initially performs a frequency 
spectrum analysis on the total data record and plots (via CALCOMP) this spectrum (Shot 
Spectra). See Volume I of this report series for an example. A computer listing provides 
the maximum spectral level and the frequency at which it occurs. The second plot presented 
4 (Shot Ampli) is the time domain data record employed in generating the initial frequency 
| spectrum. The data record is then divided into a series of small preselected time intervals 
q of equal duration and frequency spectra for each interval is computed. The resulting plot 
(Interval) is a series of frequency spectra for each time interval presented in a water-fall 
format. 


The plots length (time/frequency axis) can be varied in length from 2 inches to 28 

4 inches. All spectral data is normalized to the peak level of the initial spectral analysis. 

a The height and separation of the individual spectral plots may be varied. Highpass filtering 
: parameters (Time Domain Filter), if filtering of the data record is performed by this pro- 
gram, are also indicated on the graph. 


Ill. OUTPUT FORMATS 


The above described programs have four possible output formats. These formats 
are: 


1) Seismic Data Acquisition System (SDAS) Formatted Tape 
2) Seismic Data Display System (SDDS) Formatted Tape 
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3) Plot Tape (CalComp Tape) 
4) Data Listing 


A. SDAS FORMATTED TAPES 


This format of data output is normally used when processing the data with the follow- j 
ing programs: MAINO1, MAINO2, MAINO3, and MAINO4. This output format is used to save 
the processed data in order that other programs may process the data further. 


The SDAS formatted tape is both an input and output format. Tapes generated by the hy 

SDAS hardware are 7 track 556 Bits Per Inch. SDAS tapes generated by these programs are 
7 track 800 BPI. The first 108 bits of information, 3-36 bit comsuter words, are the | 
Header information coded in'' Binary Coded Decimal''. The program automatically decodes 
this information (see Figure 1). 
i 


B. SDDS FORMATTED TAPES 


This is an output formatted tape with 36 db dynamic range. The output format, when 
played back on the SDDS (Seismic Data Display System), will produce a 36 db dynamic range 
analog representation of the digital signal. 


The format of the tape starts with an EOF mark. Then the digital data is converted 
to a binary number. The resulting binary number is between 111 111 (all one's) and 000 000 
(all zero's). There is a'' Record Gap" after 1980 data samples (this is always the number of 
data samples per record). This format requires the number of data samples converted be 
limited to 1980 data points by: 


1) Converting the first 1980 data points (even if less data points exist) 
2) Resampling the data, reducing the information to 2000 data points, and con- 
verting the first 1980 points of the resulting data (Davis, 1970) 
Cc. PLOT TAPE (CALCOMP TAPE) 
This type of output is used to obtain a wiggle trace display of the data. The output 
tape is played back by a CalComp plotter. 


This format has a maximum physical length of 30' for each record of data plotted, 
although shorter lengths may be used. The plotter's minimum step change between each 
data point is . 005 in. 


D. DATA LISTING 


This format prints on a computer listing all data values for each record after 
processing. The data points are read sequentially from left to right. 
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Iv. PROGRAM OPERATION (UNIVAC 1108 COMPUTER) 


A. MAINO1 


The following procedure is recommended for operating the Time Domain Filter. 
First the "General Format’ is shown, then a specific sample is given. For meaning of the 
terms in the "General Format, check Glossary of Program Terms. All formats are free 


field formats, FORMAT ( ), except as noted. 


GENERAL FORMAT 


@ASG,T IN., 8C, Program Tape # 

@ CAT, P W, F40//POS/5 

@ COPY,G IN.,W. 

@FREE IN. 

@ASG,T IUNIT.,8C, DATA TAPE # 
@ASG,TM NUNIT., 8C, SAVE02 
@ASG,T 9.,8C,SAVE02 

@ASG,T (UNIT +1). ,8C, SAVE02 

@ REWIND IUNIT. 

@ XQT W.AMAINO1 

IMOVE, INCRE, NUNIT , IPRT , NDATAP 
KPASS, KOUT, 12000, AMAX, IFILE, IPOINT 
IUNIT , NUNIT 

SCALEX, SCALEY 

CUT (1), H(1), N(1) 

CUT (2), H(2), N(2) 

@TPNO NUNIT. 

@TPNO 9. 

@TPNO (UNIT +1). 


COMMENT 


NECESSARY ONLY ONCE 
PER DAY 


INPUT DATA TAPE 

IF SDDS OUTPUT SELECTED 
IF PLOT OUTPUT SELECTED 
IF SDAS OUTPUT SELECTED 


IF PLOT OUTPUT SELECTED 
IF KPASS = 2,3 
IF KPASS = 1,2 
IF SDAS OUTPUT SELECTED 
IF PLOT OUTPUT SELECTED 
IF SDAS OUTPUT SELECTED 


| SPECIFIC EXAMPLE (MAINO1) 
{| 
{| 
} @ASG,T IN., 8C, 36790 
() @CAT, P W., F40//POS 5 
hua @COPY,G IN. , 8C, SAVE02 
| @FREE IN. 
WJ : @ASG,T 10., 8C, 3505 
( @ASG,T 9.,8C,SAVE02 


@ASG,T 11., 8C, SAVE02 
@XQT W.AMAINO1 


: 1, 2,1200, 100, 4020 
. | 3,6, 2000, 2.0, 0,1 
oo 10,2 
f| 0.5,14.0 
| ms 4. 0,17. 0, 80 
| | { 40.0, 17. 0, 80 
3 @TPNO 9. 
|| @TPNO 11. 
{| 
tJ 
| 0 
im 
1 
fi 9 


B. MAINO02 


The following procedure is recommended for operating the Frequency Domain Filter. 
First the "General Format' is shown, then a specific sample is shown. For meaning of the 
terms in the "General Format’, check Glossary of Program Terms, All formats are free 
field formats, FORMAT ( ), except as noted. 


GENERAL FORMAT COMMENT | om 


@ASG,T IN. , 8C, PROGRAM TAPE # 


i @ CAT, P W., F40//POS/5 NECESSARY ONLY ONCE 
q @ COPY, GIN.,W. PER DAY 
@ FREE IN. 

@ASG,T IUNIT.,8C, DATA TAPE # INPUT DATA TAPE 
@ ASG,TM NUNIT. , 8C, SAVE02 IF SDDS OUTPUT SELECTED 
j @ASG,T 9.,8C, SAVE02 IF PLOT OUTPUT SELECTED 
@ASG,T (UNIT +1). ,8C, SAVEO2 IF SDAS OUTPUT SELECTED 


@ REWIND IUNIT. 

@ XQT W.AMAINO2 

IMOVE, INCRE, NUNIT, IPRT, NDATAP 
KPASS, KOUT , 12000, AMAX, IFILE, IPOINT 


IUNIT , NUNIT 

FRQLOW, FRQHI, KFIL,NZERO,PHASE,TSHIFT = IF BANDPASS FILTER SELECTED 
(KPASS = 1) 

@TPNO NUNIT. IF SDDS OUTPUT SELECTED 

@TPNO 9. IF PLOT OUTPUT SELECTED 


@ TPNO (UNIT +1). IF SDAS OUTPUT SELECTED 


| 
P| 
SCALEX, SCALEY IF PLOT OUTPUT SELECTED 
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| 

t 

| ; SPECIFIC EXAMPLE (MAINO2) 

j 

U 

U 

\ | @ASG,T IN., 8C, 36790 

bg @CAT, P W., F40//POS/5 i 

| (| @COPY,G IN.,W. i 

fey @ FREE IN. : 

| | @ASG,T 10., 8C, 3505 ; 

@ASG,T 9.,8C,SAVE02 fi 

| (| @ASG,T 1. ,8C, SAVE02 i 

@ REWIND OUT. q 

| {| @ XQT W.AMAINO1 i 

oe 1, 1, 1200, 100, 4010 | 

| { 1,3, 2000, 2.0, 0,1 | 

10,1 y 

i) 0.5, 14. 0 if 

Re . 100, 270, 1,0, 0,0, 0.0 ii 
i | @ TPNO 9. 
@TPNO 1. 
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Cc. MAINO3 : 

The following procedure is recommended for operating the Record Stacking Program. 
i First the "General Format" is shown, then a specific sample is shown. For meaning of the 
terms in the ''General Format, check Glossary of Program Terms. All formats are free 
field formats, FORMAT ( ), except as noted. 


GENERAL FORMAT COMMENT 


@ASG,T IN.,8C. PROGRAM TAPE # 


@CAT,P W., F40//POS/5 NECESSARY ONLY ONCE | ( 
@ COPY,G IN.,W. PER DAY a 
@ FREE IN. 

@ASG,T IUNIT.,8C, DATA TAPE # INPUT DATA TAPE E 
‘@ASG,TM NUNIT. , 8C, SAVE02 IF SDDS OUTPUT SELECTED | 
@ASG,T 9.,8C, SAVE02 IF PLOT OUTPUT SELECTED 

@ ASG,T (UNIT+1). ,8C, SAVE02 IF SDAS OUTPUT SELECTED 

@ REWIND IUNIT. 


@ XQT W.AMAINO1 
IMOVE, INCRE, NUNIT, IPRT, NDATAP 
IRRAY, KOUT , 12000, AMAX, IFILE, IPOINT 


IUNIT , NUNIT 
SCALEX, SCALEY IF PLOT OUTPUT SELECTED f 
@TPNO NUNIT IF SDDS OUTPUT SELECTED 
@TPNO 9. IF PLOT OUTPUT SELECTED | 
@ TPNO (UNIT +1). IF SDAS OUTPUT SELECTED 
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SPECIFIC EXAMPLE (MAINO3) 


@AST,T IN., 8C, 36790 
@CAT,P W., F40/POS/5 
@ COPY,G IN. ,8C, SAVE02 
@ FREE IN. 

@ASG,T 10. , 8C, 3505 
@ASG,T 9.,8C, SAVE02 
@ASG,T 11., 8C, SAVE02 
@ REWIND OUT. 

@ XQT W. AMAINO1 

1, 2, 1200, 100, 4010 

10, 6, 2000, 2.0,0,1 

10,1 

0.5,14.0 

@ TPNO 9. 

@ TPNO 11. 


D. MAIN04 


The following procedure is recommended for operating the Wiener Filter Program. 
First the General Format’ is shown, then a specific example is shown. For meaning of 
terms in the ''General Format', check Glossary of Program Terms. All formats are free 
field formats, FORMAT ( ), except as noted. 


ar 


GENERAL FORMAT COMMENT 


@ASG,T IN.,8C, PROGRAM TAPE # 


@ CAT, P W., F40//POS/9 NECESSARY ONLY ONCE 
@ COPY,G IN.,W. PER DAY 
@ FREE IN. 


@ASG,T IUNIT.,8C, DATA TAPE #1 
@ASG,T KUNIT.,8C,SOURCE SIGNATURE TAPE IF SEPARATE TAPE FOR BASIC 


@ ASG,TM NUNIT. , 8C, SAVE02 IF SDDS OUTPUT SELECTED 
@ASG,T 9.,8C, SAVEC2 IF PLOT OUTPUT SELECTED 
@ASG,T (UNIT +1). ,8C, SAVE02 IF SDAS OUTPUT SELECTED 


@ XQT W.AMAINO4 

IMOVE, INCRE, NFILE, IPRT , NDATAP 

IUNIT,, NUNIT 

KPASS, KOUT , IEXPND, 12000, AMAX, FILE, IPOINT 
ITER, LENGTH, INCLEN, LS, LF 

IAVE, KMOVE, KUNIT 


SCALEX, SCALEY IF PLOT OUTPUT SELECTED 
LA, LB, LD, IST1, KB, KE 

D@) FORMAT 10F 8.5 

| IFTD IF KPASS=1, 2,3 

CUT (1), H(1), NQ) IF KPASS=2, 3 

CUT (2), H(2), N(2) IF KPASS=1, 2 

| FRQLOW, FRQHI, KFIL,NZERO, PHASE,TSHIFT _IF KPASS-5,6 

| BASIC WAVESHAPE (CARD FORMAT 10F8. 5) IF IEXPND IS -1 


FRQLOW, FRQHI, KFIL, NZERO, PHASE, TSHIFT IF KPASS=6, 7 
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U | 
: l SPECIFIC EXAMPLE (MAIN04) : 
{| @ASG,T 10., 8C, 7651 
@ASG,T 20., 8C, X0001 

| (| @ASG,T 9. ,8C, SAVE02 | 
| _ 1,1, 1000, 100, 4010 
 { 10,1 
es 4,14, 0,0,2.0,0,1 ' 

| 1, 160, 6, 102, 140 

| U 30, 130, 20 
0.5,14.0 : 

Ms 240, 160, 100, 1, 400, 400 ; 

(| DESIRED WAVESHAPE (CARD FORMAT 10F8.5) | 


IM, AS eV PE 


j 
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E. MAINO7 


The following procedure is recommended for operating the Desired Waveshape 
Program. First the General Format is shown, then a specific example is shown. The 
meaning of the terme in the General Format" is given in the Glossary of Program Terms. 
All formats are free field formats, FORMAT ( ), except as noted. 


GENERAL FORMAT COMMENT 


@ ASG,T IN. ,8C, PROGRAM TAPE # 

@ CAT,P W., F40//POS/9 

@COPY,G IN.,W 

@ FREE IN. 

@ XQT, W. AMAINO7 

FSAMP, NP, PHASE, NZERO, FRQLOW, FRQHI, KFIL, IPUNCH 

ZERO1 (1), ZERO2(1) IF ZERO IS GREATER THAN ONE. 
ONE FOR EACH FREQUENCY TO 
BE ZEROED 
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SOT, FEC, HI 
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1. 000, 1024, 0.0, 0,.010,.170, 2,0 


@ XQT W.AMAINO7 


& 
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F. SDADS 


The following procedure is recommended for operation of the Spectrum Analysis 
Program. First the General Format" is given, then a specific example. For the meaning 
of terms in the "General Format’, refer to the Glossary of Terms. All formats are free 
field, FORMAT ( ), except as noted. 


GENERAL FORMAT COMMENT 


@ ASG,T IN., 8C, PROGRAM TAPE # 

@ CAT,P W., F40//POS/9 

@ COPY,G IN.,W. 

@ FREE IN. 

@ ASG,T IUNIT,8C, DATA TAPE # 

@ ASG,T 9.,8C, SAVE02 

@ XQT W.AMAINI5 

IGO, XFAC, IFIL, ITD, ITEST , IVER 

SAMP1, SAMP2, SAMP3, TIME, ISHIFT , SHIFT 
L1, L2, DSH, IFFT 

ROF1, IFIRST, PD, KB, KE 

CUT (1), H(1), N(1) IF IFIL EQUAL 1 (HIGHPASS) 
IMOVE, INCRE, NFILE, IPRT, SCALE, IUNIT rahe. 

SMSEC, NSAMP, ITER, SCTSTK, DELBSP, DBRANG, TAPE 


IF DATA IS READ FROM SDAS 
TAPE 


DATA POINTS (CARD FORMAT 10F8. 4) IF ITEST GREATER THAN 1 
FOR ANALYSIS (TEST DATA POINTS) 
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| | SPECIFIC EXAMPLE (SDADS) 
iT 
t 
: I @ASG,T 10., 8C, 7651 
@ ASG,T 9.,8C, SAVE02 
| ] 0, 0.0,0,0,0,0 
: 4.0,1.0,1.0,2.0,1,.001 
I 0,0, 0,0 
0,1,0,1, 2000 
| | 1,1,10,1,1.0,10 
0.0,0,1,1,0,60.0,1.0 
if 
[ 
| 
{ 
{ 
if 
| E 
in! 
iF 
: 
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V. PROGRAM OPERATION (CDC 6600) 


The seismic processing programs described in this report have been designed for 
use primarily on a Univac 1108 computer. To operate these programs on the CDC 6600 com- 
puter the statements controlling input/output of data must be changed to their equivalent 
CDC 6600 statement form. These statements are shown below. 


Univac 1108 Statement Form CDC 6600 Equivalent Statement Form 
NTRAN BUFFER 
FOLD (used to unpack computer words) UNPAKF 
FOLD (used to pack computer words) ENCODE 
READ(5,190) (Unformatted READ) READ (5, *) 


The above statements are not one for one substitutions and the programs must be changed 
accordingly to accept them. The changes to be made to MAINO] are presented. Similar 
changes must be made to the other main programs. All of the subroutines that must be 
changed for operation of any of the programs are also shown. The subroutines that were 
changed are SUB001, SUB002, SUB003, SUB044. No other subroutines need be modified for 
the conversion to CDC 6600. The method of operating MAINO1 on CDC 6600 is also shown. 
The operation of the other programs is similar. 


Refer to Appendix A listing for line number reference of UNIVAC statement. 
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7 7 23 . I at ed ied renee ak 
‘s 
ra " roe = . " Diana . . " 


96 


100 
158 
159 


CHANGES TO MAINO1 


Univac Statement 


READ (5, 190) 

READ (5,190) 

READ (5, 190) 

CALL NTRAN (UNIT, 8, IMOVE) 
IF (PLOT) READ (5,190) 

IF (LOW) READ (5,190) 

IF (HIGH) READ (5,190) 

IF (SDDS) CALL NTRAN (NUNIT, 9) 


IF (SDAS) CALL NTRAN 
(IUNITT 1, 9) 


EXIT = 0 
IF (SDDS) CALL NTRAN (NUNIT, 9) 
IF (SDAS) CALL NTRAN (UNIT1, 9) 


ed aa ed 


122 


123 


124 
125 
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CDC Statement 


PROGRAM (INPUT, OUTPUT , TAPE5= 
INPUT, TAPE6=OUT PUT, 
TAPE1, TAPE9, TAPE10, 
TAPE11) 


READ (5, *) 

READ (5, *) 

READ (5, *) 

CONTINUE 

IF (PLOT) READ (5, *) 

IF (LOW) READ (5, *) 

IF (HIGH) READ 6, *) 

IF (. NOT.SDDS) GO TO 123 


BUFFER OUT (NUNIT,1) (DATA(1), 
DAT A(2)) 


IF (UNIT (NUNIT)) 122, 122, 122 

REWIND NUNIT 

END FILE NUNIT 

IF (. NOT.SDAS) GO TO 125 

BUFFER OUT (IUNIT1,1) (ATA(1), 
DATA(2)) 

IF (UNIT (IUNIT1)) 124, 124, 124 

END FILE IUNIT1 

CONTINUE 

EXIT =0 

IF (SDDS) END FILE NUNIT 

IF (SDAS) END FILE IUNIT1 

KK = NDATAP 


B. CHANGES TO SUB001 


Line Univac Statement CDC Statement 


. 14 DIMENSION INDATA (32,2), DIMENSION INDATA (1200), IST (2), 
i IST (2), IHEAD (3), ITHEAD(3), IDATA (27) 
IDATA G7) EQUIVALENCE (IHEAD@), INDATA(1)) 
20 CALL NTRAN (NUNIT, 9) ENDFILE NUNIT 
j 23 CALL NTRAN (UNIT, 2,3, IHEAD IF (KK. GT.6000) KK=6000 
Oh RE) NK = KK/5+2 
| BUFFER IN (IUNIT,1) (IHEAD(1), 
i INDATA(NK)) 
| 24 IF (IST(K).EQ.-1) CALL 
NTRAN (UNIT, 22) IF (UNIT(UNIT)) 40, 31, 32 
25 IF (IST(K).GT.0) GO TO 40 CONTINUE 
26 31 KK =IST(K) KK = -2 
GO TO 33 
32 KK=-3 


33. CONTINUE 


REMOVE CARDS 37 THRU 55 REPLACE 
THESE CARDS 
CXXXX KK IS THE NUMBER OF SAMPLES 
READ 


KK =0 | 
CXXXX READ DATA RECORD = 
44 CALL LENGTH X(UNITT, IR, NUBC) 
IR = IR-2 
NR = NUBC/12 
KK = IR*5 - NR 
IF (PR.EQ.1) PRINT 862, IR, NUBC, “4 
NR, KK 


a 862 FORMAT (4110) 


7 CALL UNPACK (HEAD(2), 12,5, 4 
a JDATA) - 
b. IF (JDATA(5). GE. 4000B) JDATA(5)= F 
; 4000 - JDATAG) 4 
a 


q 22 Ny 
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CDC Statement 

DATA (KK) = JDATA(5) 
NK = NK-2 

DO 63 I =1, NK 


CALL UNPAKF (INDATA()), 12,5, 
JDATA) 


DO 64 J =1,5 


IF (JDATA(J). GE. 4000B) JDATA(J) = 
4000B - JDATA(J) 


KK = KK+1 
DATA (KK) = JDATA (J) x . 0048852 

64 CONTINUE 

63 CONTINUE 


i 
a 
a 
4 
i 
| 
I 
| 
| 
{ 
I 
I 
| 
[ 
| 
[ 
[ 
§ 
qi 


Cc. CHANGES TO SUB002 


Univac Statement 


17. DIMENSION DATA(1), NOUT(330) 


23 DO 20J =1,330 

25 DO10K=1,6 

34 FLD (IF1,6,NOUT(J)) = FLD 
(30, 6, NDATA) 


35 IF1=IF1+6 


47 DO 35 J = 1,330 

49 DO 30 IX = 1,6 

58 FLD (IF1,6, NOUT (J)) = FLD 
(30,6, NDATA) 


59 IF1=IF1+6 


62 38 CALL NTRAN (NUNIT, 1,330, 


NOUT (1), ISTAT , 22) 


63 IF (ISTAT.EQ.-1) CALL NTRAN 


(NUNIT, 22) 


64 IF (STAT.LT.1) GO TO 40 
71 RETURN 


| D. CHANGES TO SUB003 


Untvac Statement 
11 DO 411#1,3 
7 (REMOVE 12 - 18) 
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39 


41 


CDC Statement 

DIMENSION DATA(1), NOUT (208) 
DO 20 J =1,198 

DO 10 Xx = 1,10 

Ty = 198 + 1X 

NOUT (TY) = NDATA 

CONTINUE 


CALL ENCODE (10, 1, NOUT(J)) 
NOUT (199) 


FORMAT (10R1) 
DO 35 J = 1,198 
DO 30 IX = 1, 10 

TY = 198 + 1X 

NOUT (TY) = NDATA 
CONTINUE 


CALL ENCODE (10, 2, NOUT(J)) 
NOUT (199) 


BUFFER OUT (NUNIT, 1) (NOUT(1), 
NOUT (198) 


ISTAT = -2 

IF (UNIT (NUNIT )) 39, 40, 41 
CONTINUE 

RETURN 

ISTAT =-3 


CDC Statements 
CALL UNPAKF (IHEAD, 4, 27, IDATA) 


ee 


Sia os 


ee 


20 
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CHANGES TO SUB044 


Univac Statement 


DIMENSION DATA (1), NOUT (332) 
IHEAD@) 


CALL NTRAN (NUNIT, 1,3, HEAD 


(1), ISTAT , 22) 


IF (ISTAT.EQ.-1) CALL NTRAN 
(NUNIT , 22) 


IF (ISTAT.LT.1) GO TO 40 
JI=0 


IBUF = 332 
ITEST = 332 


REMOVE CARDS 22 thru 25 


28 
34 


35 
36 
37 


39 


41 


45 


47 


DO 20 IX =1,3 
IF (DATA(JI). LT. 0) ISIGN = 1 


FLD (IF 1,1, NOUT(J)) = ISIGN 
IF1=IF1+1 


FLD (IF1,11, NOUT(J)) = FLD 
(25,11, NDATA) 


20 IF 1=IF1+11 


CALL NTRAN (NUNIT, 1, IBUF, 
NOUT (1), ISTAT, 22) 


IF (ISTAT.EQ. -1) CALL NTRAN 
(NUNIT , 22) 


IF (ISTAT.LT.1) GO TO 40 
IF (KOUNT.LT.1) RETURN 
ITEST =ITEST + 332 

GO TO 10 


25 


20 


39 
41 


CDC Statement 


DIMENSION DATA(1), NOUT (1207), 
IHEAD(3) 


J = IHEAD(2) 


CALL ENCODE (10, 1, IHEAD(2))I, 
IDATA(1) 


FORMAT (A8, R2) 

JI =2 

NOUT (1) = IHEAD(1) 

NOUT (2) = IHEAD(2) 

IBUF = KK/5 

IF (KK.GT.6000) KK = 1200 


DO 20 Ix =1,5 


IF (DATA(JI). LT.0) NDATA = 4000B 
+ NDATA 


Ty = 1200 + Ix 
NOUT (IY) = NDATA 
CONTINUE 


CONTINUE 
ENCODE (10, 2, NOUT(J)) NOUT (1203) 
FORMAT (5R2) 

IBUF = IBUF +2 


BUFFER OUT (NUNIT, 1) (NOUT(1)), 
NOUT (IBUF) 


ISTAT = -2 


IF (UNIT (NUNIT)) 39, 40, 41 
RETURN 

ISTAT =-3 

CONT INUE 
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F. | MAINO1 PROCEDURE | 


The following procedure is recommended for operation of the Time Domain Filter. 
First the ''General Format" is given, then a specific example. For the meaning of terms in 
the ''General Format", refer to the Glossary of Terms. 


GENERAL FORMAT COMMENT 
X862, F4, MTa. B. E. ECKSTEIN a =1 IF NO SDAS OUTPUT 
ATTACH LIBRARY, ID=NB, SN=USET1. a =2 IF SDAS OUTPUT ry 
UPDATE Q, P=LIBRARY. od 
FTN, I, B=BEE, SL. 
REQUEST, TAPE1, *PF,SN=USET1. IF SDDS OUTPUT SELECTED 
LABEL, TAPE10, R, L=SDAS REEL #1615 INPUT DATA TAPE 
LABEL, TAPE11, W, L=SDAS. , RING, X=SV, T=b. IF SDAS OUTPUT SELECTED 
SKIPF TAPE10. 


ROUTE TAPE 9, DC=PU, TID=AF, DEF, EC=80COL. IF PLOT OUTPUT SELECTED 


ATTACH, PLOT, ZETAPLOT , ID=MZ, SN=USZT2, 
MR=1. IF PLOT OUTPUT SELECTED 


LDSET PRESET = ZERO. 

LOAD BEE. 

LOAD PLOT. IF PLOT OUTPUT SELECTED 

EXECUTE. 

CATALOG TAPE1 IF SDDS OUTPUT SELECTED 7] 
| ‘ a 
7 *C MAINO1 
: ? 

IMOVE, INCRE, NUNIT, IPRT , NDATAP 

KPASS, KOUT , 12000, AMAX, IFILE, IPOINT 


PIE OTT Pye ET, OS PET ITE IO) 4 Re 
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TUNIT , NUNIT 8, 
: } 
3 SCALEX, SCALEY IF PLOT OUTPUT SELECTED 

: CUT (1), H(1), NQ) IF KPASS = 2,3 , 


CUT (2), H(2), N(2) IF KPASS = 1,2 


OR TERT A eT PNT 


aq an cf & a & & 


SPECIFIC EXAMPLE (MAIN0O1) 


x862, F4, MT2. B. E. ECKSTEIN 
ATTACH LIBRARY, ID=NB, SN=USET1. 

UPDATE Q, P=LIBRARY. 

FIN, I, B=BEE, SL. 

LABEL, TAPE10,R, L=SDAS. 

LABEL, TAPE11, W, L=SDAS. , RING, X=SV, T2364. 
SKIPF TAPE10. 

LDSET PRESET = ZERO. 

LOAD BEE. 

EXECUTE. 

? 

*C MAINO1 

? 

1, 2, 1200, 100, 4000 

3, 8, 2000, 2.0,0,1 

10,1 

4,0,17. 0, 80 

40.0, 17.0, 80 

? 

? 
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VI. SUMMARY 


Six software programs have been developed to enhance seismic data quality. The 
Time Domain Filter is a zero phase Linnette Filter, The Frequency Domain Filter filters 
seismic data using rectangular, cosine, gaussian or triangular shading of the frequency con- 
tent. The vertical record stacking program sums multiple records having the same source 
to hydrophone array distance, but different common depth points. The Wiener Filter pro- 
gram minimizes the power existing between a desired output and an actual output signal re- 
sulting in the least squares optimum filtering of the seismic data. The desired waveshape 
program generates a zero phase band limited pulse with predetermined spectral frequency 
characteristics. 


The spectral analysis program analyzes the frequency content of a time series. The 
output formats of these programs allow either display of information or further processing 
of data. The SDAS Formatted tape output, a 72 dB dynamic range data representation, allows 
further processing of the data to occur before display. The SDDS Formatted tapes, a 36 dB 
dynamic range data representation, is converted to an analog signal by the SDDS. The Plot 
Tape contains wiggle trace display formatted data. The computer data listing may contain 
a list of each data point in each record. 


At present, these programs can operate on the UNIVAC 1108 or the CDC 6600. 
Because of the amount of computer time required by operation of these programs, they 
should be considered for use on a high speed computer such as the TI-ASC. 


VII. GLOSSARY OF PROGRAM TERMS 


Term Meaning 
AMAX Determines the normalizing value of the data for SDDS and PLOT outputs. 


If AMAX is less than 0, the program assumes the maximum value of the 
first record for the normalizing value and uses the updates given by SUB002; 
if AMAX equals 0, the program assumes the maximum value of the first 
record is the normalizing value; if AMAX is greater than 0, the program 
uses the value of AMAX as the normalizing value. 


BY) Is the array in which the transmitted pulse waveform is stored for use by 
the Wiener Filter Program. The waveform is called the BASIC WAVESHAPE 
and has LB number of points. This is card input data if IEXPND is less 
than 0. (Card format is 10F8.5. ) 


CUT(D) Is one of the input parameters for the Time Domain Filter. The value re- 
lates to the cutoff frequency but the actual cutoff frequency is CUT (1) plus 
H() (Value in HZ). I=1 for lowpass parameters, I=2 for highpass parameters. 


Dq@) Is the array in which the DESIRED WAVEFORM, used by the Wiener Fil- 
ter, is stored. The waveform normally is developed by MAINO7. (Card 
format is 10F8.5. ) 


DBRANGE Determines the range in db of the graph output of the Spectral Analysis 
Program, 
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DDATA(]) Is the array input data stored by the Spectral Analysis Program. ITEST 
number of points must be supplied to the program if analysis is to be per- 
formed on other than SDAS tapes. 


DELBSP Is used by the Spectral Analysis Program when picking the first return. It 5 
moves the first return point (chosen on peak value) DELBSP points earlier. he 
This value should be zero if program is not choosing the first return, 


FRQHI The highcut filter frequency used in the FREQUENCY FILTER and the e 
DESIRED WAVESHAPE programs (Value is in KHZ). i 


FRQLOW Is the lowcut filter frequency used in the FREQUENCY FILTER and the 
DESIRED WAVESHAPE programs (Value is in KHZ). 


SS /_- — = fe 


FSAMP This is the sample rate of the DESIRED WAVEFORM (Value is in KHZ). 


[ H(@) The slope parameter for the Time Domain Filter (Zero is maximum slope. ) 
[ I=1 for lowpass parameters, I=2 for highpass parameters. 

: IARRAY Is the number of records to be stacked. 

i IAVE Is the number of records stacked (added together), by the Wiener Filter 


developing the BASIC WAVESHAPE. This does not affect the wiener filtering 
of data records with filter weights directly. 


STAY BNO RNY IAAT ING ENT AEN NT NONE SN Ib OAL nada 


Becta d 
> ¢ 


IDSH If not equal to 0, causes printout of the first IDSH values moved by the sub- 
routine DAT SHI. 


* 
a 
’ 

> 


[ 
I 
[ 
l 
! 
i 
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IEXPND Causes the Wiener Filter program to change processing of the BASIC 
WAVESHAPE. If IEXPND is less than zero, read in the basic waveform; 
if IEXPND equals 0, do not expand the bubble period; if IEXPND is greater 
than 0, expand the bubble period according to the input parameter length. 


IFFT If not equal to zero, causes printout of the first IFFT values of output by the 
subroutine FFT (the subroutine which converts data from the time domain 
to the frequency domain). 


IFIL Is used by the Spectrum Analysis, if equal to 1, causes highpass filtering 
of the data by the Time Domain Filter. It also demands the reading of the 
Time Domain Filter parameters. 


IFILE Changes the gain of the digital data for the first IFILE number of records. 
If IFILE is positive, it multiplies the data by 3; if IFILE is negative, it 
divides the data by 3. 


IFTD Is used in MAINO4 when processing data with the Time Domain Filter. If 
IFTD equals 1, filters input SDAS data used before creating filter weights; 
if equal 2, filter all input SDAS data; if equal 3, filters input SDAS data used 
to convolve with filter weights. 
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IMOVE 


INCLEN 


ISHIFT 


ITER 


ITER1 


ITEST 


IUNIT 


IVER 


12000 


KB 


Is the first file processed by the program. 


Is used in MAIN04 when repetitive expansion of the bubble pulse period is 
desired. The bubble pulse period is increased by INCLEN data samples 
each iteration. 


Indicates which records are read (every increth record). This value must 
be greater than zero. 


Indicates output starts at the [POINT point in the data array. 


Indicates how often Header Records and PLOT Indication is printed on list- 
ing. The first Header Record and PLOT Indicator is always printed. 


Punch output on card. 


Used by Spectrum Analysis Program, shifts the data in time by this number 
of data samples. 


Shifts the desired waveshape by IST1 data points. This value must be 
greater than -1 but less than 30 if LD equals 100. 


Indicates that the interval spectral analysis will be bandpassed analyses 
plotted in the time domain. 


Indicates the building of filter weights should be repeated ITER number of 
times. Used for testing with bubble pulse expansion, 


Is used in the Spectrum Analysis Program. When ITD equals 1, Low passes 
interval analysis data rather than bandpass data. 


Is a control of Spectrum Analysis Program. If ITEST equals 1, one point is 
entered at the 300th data point, sample rate is 1000, sample length equals 

1 sec, and analysis is performed; if ITEST is greater than 1, ITEST values 
are read in, sample rate set to 1000 and sample length is set appropriately 
and analysis i¢ performed; if ITEST equals 0, analyzed data is read from 
SDAS tape. 


Is the file name of the input SDAS tape unit. This is a number not a name! 


Indicates that the Interval Analysis by Spectral Analysis Program is plotted 
vertically if IVER equals 1. 


If equals 1, the data is limited to 2000 data points; {f 12000 is greater than 1, 
the data is limited to 12000 data points; if equal to 0, uses actual number of 
data points. 


Is the first data point observed when determining the maximum value of 
the data. 


Is the last data point observed when determining the maximum value of 
the data. 

Is the type of taper applied by the Frequency Domain Filter. 

KFIL =0 ino taper 

KFIL =1 cosine taper 

KFIL=2 Gaussian taper 

KFIL=3 Triangular taper 


Is the record used by MAIN04 containing the BASIC WAVESHAPE. 


Determines the output type. 
KOUT= SDAS SDDS 


NO NO 
YES YES 
YES YES 
YES YES 
YES 
NO 
NO 
NO 
NO 


0 
1 
2 
3 
4 
5 
6 
7 
8 
9 


Determines type of filtering 
If used in Time Domain Filter 


KPASS = 1 lowpass filter 
KPASS = 2 highpass filter 
KPASS = 3 bandpass filter 
KPASS = 4 no filtering 


If used in the Frequency Domain Filter 


KPASS = 1 filter 
KPASS = 0 no filter 


If used tn MAINO4 


KFIL =1 Lowpass Filter Time Domain 
KFIL = 2 Highpass Filter Time Domain 
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KFIL =3 Bandpass Filter Time Domain ; 
KFIL = 4 No filtering 


KFIL =5 Frequency Domain Filter section 1 
KFIL =6 Frequency Domain Filter both sections 
KFIL =7 Frequency Domain Filter section 2 
KUNIT Is BASIC WAVESHAPE tape unit file number. 
Li 
LA Is the number of wiener filter weight points. 
LB Is the number of BASIC WAVESHAPE points (before bubble pulse expansion). 7 
LD Is the number of DESIRED WAVESHAPE points (before shifting by IST1). ie 
\4 
LENGTH Is the initial length of the bubble pulse period after expansion of the BASIC L 
WAVESHAPE, 
4 
LF Determines the last data point of the BASIC WAVESHAPE for bubble pulse 


expansion (relative to first data point read from tape). 


LS Determines the first data point of the BASIC WAVESHAPE for bubble pulse 
expansion (relative to first data point read from tape). 


Ll Controls the plotting of the line on the computer listing if 0, do not list plot; 
if greater than 0, list L1 points in plot (this does not affect plot output). 


L2 Controls the plotting of the first return on computer listing. If 0, do not 
list plot; if greater than 0, list L1 points starting at picked first return. 


N(), NIA These two are the same parameter. They control how many filter weights 
are used by Time Domain Filter for filtering. N(I) equals 1/2 filter weights. 
I=1 for lowpass parameters, I=2 for highpass parameters. 


NDATAP The maximum number of data points a data tape record will have. It may be 

greater than the actual number of data points. 
: NFILE The last record read. 

NSAMP The number of "interval spectral analysis" performed by the Spectral 
Analysis Program. 

NUNIT This is the output unit for SDDS output tape and normally is unit 1. Tape 
must be assigned medium density. 

NZERO Indicates the number of rejection filters used by the Frequency Domain | { 
Filter. 


a PHASE Shifts the phase of the data being filtered. Normally this value is zero. 


SAMP1 Determines the height of the graphs plotted for the frequency analysis of all 
data and the time domain plot of the data. 

SAMP2 Determines the offset between successive interval analyses. If SAMP2 
equals 0, then SAMP1 = SAMP2 = 10/NSAMP. 

SAMP3 Determines the height of interval analysis. 

SCALE If analog tapes are digitized, but at a different than real time rate, this 


factor compensates for the change. SCALE = speed at playback 
speed at recording 


SCALEX Vertical height (in.) of the plot information. 

SCALEY Horizontal length (in. ) of the plot for 2000 points of information. 

SCTSTK Value determines how many records are stacked. 

SHIFT The first point plotted in the time domain spectrum. 
'SMSEC "Interval Spectral Analysis" number of data points or frequencies for 

analysis. 

TAPE This number is written on the graph which is plotted. 

TIME The last data point, relative to time, plotted in the "time domain spectrum". 
TSHIFT Time shifts the data by TSHIFT number of points. 

XFAC Scales down the size of the plot in the Spectrum Analysis Program. 
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